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How much does the curvature perturbation change after it leaves the horizon, and when 
should one evaluate the power spectrum? To answer these questions we study single field 
inflation models numerically, and compare the evolution of different curvature perturba- 
tions from horizon crossing to the end of inflation. We find that e.g. in chaotic inflation, 
the amplitude of the comoving and the curvature perturbation on uniform density hy- 
pcrsurfaces differ by up to 180 % at horizon crossing assuming the same amplitude at the 
end of inflation, and that it takes approximately 3 efolds for the curvature perturbation 
to be within 1 % of its value at the end of inflation. 
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1. Introduction 

It is well known that the curvature perturbations on both uniform density hyper- 
surfaces, C, and on comoving hyper surf aces, 1Z, are conserved on large scales. The 
standard approach used to calculate the power spectrum of perturbations after hori- 
zon crossing assumes that the limit k — > has been reached, see RefP However, 
immediately after horizon crossing the wavenumber will not yet have become suffi- 
ciently small for this limit to be accurate and gradient terms will still play a role. 
Exactly how long this evolution will last and the size of the errors if the curvature 
perturbation is evaluated too early are issues addressed in this work. 2 

2. Equations and Numerics 

We consider single field inflation in a homogeneous, isotropic, Friedmann-Robertson- 
Walker background, with scalar perturbations of the field and metric to first order, 
and work in flat gauge and Fourier space. The resulting first order Klein-Gordon 
equations are solved numerically by following Salopek et al. 3 Initial conditions for 
the background are selected for each potential following HustorP and at early times 
we assume the Bunch-Davies vacuumPEI We map the initial scalar field fluctuations, 
onto conserved quantities, that remain constant in the limit k — > for adiabatic 
perturbations. We focus on the curvature perturbation on uniform density hyper- 
surfaces £ = ip + -j^8p, and the comoving curvature perturbation TZ = ip+ -^rS(p (see 
e.g. Ref.[5]). These two gauge-invariant curvature perturbations, are related by a con- 
straint equation and ( + 1Z becomes small on super- horizon scales. We show plots for 
the modes ki = 2.77 x fO~ 5 Mpc _1 , k 2 = 2.00 x IQ^Mpc" 1 (WMAP pivot scale) 
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and /v3 = 1.45 x 10 -1 Mpc -1 . The results below are for the potential U = ^m 2 (p 2 
with m = 6.32 x 10 _6 MpL but we also obtain results for U — Uo + ^m 2 if 2 , U = jXtp 4 
and U — aip 2 / 3 . For full details see Nalson et alPOur numerical results have been 
verified with a second numerical program, pyflationP 

To deduce an analytic expression for the curvature perturbation, we follow Ref. [7] 
The Eq. (3.34) given there is only valid in the large scale limit, however, it must be 
evaluated exactly when the corresponding mode crosses the horizon. 

3. Results 



As expected, we can see from Fig. 1(a) that a short time after horizon crossing there 
is no longer any appreciable evolution in either the power spectrum of C, *P((k) or 
the power spectrum of TL, V-jzik), but that there is some evolution immediately after 
horizon crossing. We find that despite £ and 1Z being equivalent very far outside 
the horizon, the difference between \1Z\ and |£| at horizon crossing can be as much 
as 20% and £ remains significantly larger for at least a couple of efolds. We also 
find that the error in evaluating the power spectra numerically at horizon crossing 
rather than using the correct analytic expression or the full numerical solution at 
late times can be as much as 180% for V^(k) and 100% for V-r (this is expected 
analytically^. Lastly we showed that to evaluate the power spectra without using 
the analytic expression, one would need to wait 3.2 and 2.9 efolds to ensure the 
answer for Vn{k) and P^, respectively, are correct to within 1% of the value at the 
end of inflation. Similar results were obtained for the other potentials. 
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(a) Percentage difference between V{k) evalu- 
ated, numerically, at the end of inflation and 
at each time step, plotted for k\ (black, left), 
&2 (blue, middle) and k$ (red, right). 



(b) Numerical solutions compared to correct 
analytic solution and naive evaluation of the 
power spectra at horizon crossing, plotted for 
fc2 (left lines) and k-j, (right lines). 



Fig. 1. The evolution of V-jz(k) and (fe) are plotted against the number of efolds, N '■ 

In Fig. [T] we compare the correct analytic solution and the naive calculation 
of the power spectrum at horizon crossing with the numerical solutions. When we 
compare the numerical solution to V-fi(k) evaluated at horizon crossing we find that 
there is a 100% error in our answer. We also see that if one were to use the analytic 
expression, but to evaluate it 'some efolds after horizon crossing' rather than at 
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horizon crossing one would underestimate the amplitude of the power spectrum. 
This corresponds to following the red line in Fig. 1(b) e.g. evaluating the analytic 
expression 4 efolds after horizon crossing incurs a 10% error. Lastly we see that 
the analytic and numerical expressions do not agree with each other shortly after 
horizon crossing, one must wait at least 3.2 efolds for these two values to agree. 



4. Conclusions 

We have quantified the evolution of the curvature perturbations after inflation and 
highlighted possible errors which can occur. As we are entering an era where we 
hope to constrain parameters to within a percent using e.g. Planck data, it is of 
particular importance that these errors are both minimised and quantified. We show 
the difference between analytic and numerical expressions close to the horizon. The 
numerical results, instantaneous values of the power spectrum at horizon crossing, 
while not of observational significance, are useful in many ways, e.g. as initial con- 
ditions for other analytical or numerical schemes operating outside the horizon. If 
we are interested in late time values we should not evaluate numerical results at 
horizon crossing, as unlike the analytic results they will not be accurate. If we are 
interested in the instantaneous values at or close to horizon crossing, the analytic 
expressions are no longer valid and one must use numerical methods. For example, 
this is important if there is a second phase of evolution which starts to dominate 
during the first three efolds after horizon crossing. In conclusion, we have high- 
lighted that confusion between the different curvature perturbations, how they are 
evaluated and when each expression is valid can introduce additional errors when 
comparing theoretical results with observations. 
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